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ABSTRACT 


A  common  type  of  digital  communication  system  is  binary  frequency  shift  keying 
(FSK)  whereby  every  T  seconds  the  transmitter  sends  a  pulse  of  one  of  two  frequencies. 
The  receiver  makes  a  decision  (every  T  seconds)  as  to  which  frequency  was  trans¬ 
mitted.  A  sub-optimum  receiver  for  this  case  obtains  estimates  of  the  two  noise  wave¬ 
forms  by  passing  received  signals  through  filters  centered  at  the  sending  frequencies 
and  then  croos-correlates  these  estimates  with  the  received  waveform.  Two  slightly 
different  versions  of  this  cross-correlator  were  considered,  and  the  probability  of 
error  for  each  case  was  calculated.  The  results  seem  to  agree  with  previous  experi¬ 
mental  work  by  Cossette  and  Wolf. 
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1.  INTRODUCTION 


A  common  type  of  digital  communication  system  is  binary  frequency  shift-keying 
(FSK)  whereby  every  T  seconds  the  transmitter  sends  a  pulse  of  one  of  two  frequencies. 
The  receiver  then  makes  a  decision  (every  T  seconds)  as  to  which  frequency  was  trans¬ 
mitted.  The  method  whereby  the  receiver  makes  this  decision  for  various  types  of 
communications  channels  is  the  subject  of  this  report. 

The  simplest  model  for  a  communications  channel  assumes  that  the  input  to  the 
receiver  is  an  attenuated  version  of  the  transmitted  signal  corrupted  by  additive 
Gaussian  white  noise  (with  zero  mean  and  power  spectral  density  S(w)  =  No/2  for  all 
a)).  This  model  has  been  thoroughly  analyzed  in  the  literature^1)  and  the  receiver 
structure  which  leads  to  t.he  minimum  probability  of  error  is  known.  Specifically,  let 
us  assume  that  the  received  waveform,  r(t),  is  given  as 

\2E/T  sin  (u)ot  +  9q) 
or 

VSeTt  8in  (o^t  +  Gj) 

where 

(a)  cdoT  =  k2  it  k,  an  integer 

(b)  GUjT  =  p2rr  p,  an  integer  not  equal  to  k 

(c)  E{n(t)}  =0,  E{n(t)n(t -t)}  =  6(r) 

and 

(d)  The  a  priori  probabilities  of  each  sinusoid  are  equal. 

Many  receiver  structures  can  be  given  all  of  which  have  identical  performance. 

Two  receiver  structures  which  lead  to  the  minimum  probability  of  error  for  this 
simple  channel  model  are  given  in  Figure  1.  The  probability  of  error  for  these  re¬ 
ceivers  is: 


(2) 


1 


The  above  situation,  termed  coherent  detection,  may  not  be  realistic  for  many 
reasons.  One  of  these  reasons  is  that  it  is  assumed  that  the  receiver  has  knowledge  of 
the  exact  phase  of  both  sinusoids .  If  the  receiver  knows  nothing  of  the  phase  of  the 
received  sinusoids  and  it  is  either  undesirable  or  impractical  to  assume  that  it  has 
estimated  this  phase  from  previous  pulses,  then  we  can  add  the  additional  restriction 
that: 


(e)  The  phases  and  #2  are  independent  random  variables,  each  having  a  proba¬ 
bility  density  function  which  is  uniform  over  the  interval  (0,  2 it). 

Again  the  optimum  receiver  which  leads  to  the  minimum  probability  of  error  is 
known.  Two  such  optimum  receiver  structures  are  given  in  Figure  2.  The  probability 
of  error  for  these  receivers  is 


(3) 


Continuing  with  the  idea  of  making  tne  mathematical  model  of  the  channel  more 
general  so  that  it  applies  to  a  wide  class  of  channels,  it  is  now  assumed  that  there  are 
statistical  fluctuations  in  the  amplitude  of  the  FSK  signals.  The  simplest  situation  to 
consider  is  the  case  where  the  amplitude  \2E/T  is  constant  over  any  one  pulse  period 
but  varies  in  a  statistical  manner  from  pulse  to  pulse.  In  that  case,  the  previous  re¬ 
ceiver  structures  are  still  optimum  but  the  average  probability  of  error  is  given  as 

.00 

=  /  (probability  of  error  without  fading)  p(E)dE  (4) 

o 


where  p(E)  is  the  probability  density  function  of  the  signal  energy  E,  which  is  a  random 
variable.  Amplitude  fluctuations  are  usually  accompanied  by  an  unknown  phase  so  (4) 
becomes : 


00  _JL_ 

—  r  l  2No  1 

pe  -  J  2  e  p(E)dE  =  Mg  (  -  )  (nonchoerent  case)  (4a) 

where  Mg  (Jv)  is  the  characteristic  function  of  the  signal  energy  "E‘\ 

If,  however,  the  phase  is  known  (or  a  very  ;jood  estimate  is  made  of  the  phase)  then 
(4)  becomes 


p(E)dE  (coherent  case) 


(4b) 
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r(f) 
- ■» 


r(t  | 


Figure  2. 


(o)  QUADRATURE  -  CORRELATOR  RECEIVER 


<  b  )  MATCHED  FETE  Hr  -  ENVELOPE  DETECTOR  RECEIVER 


Two  Realizations  for  Optimum  Receiver  for  Additive  Gaussian  Noise  and 

No  Phase  Information 
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A  common  type  of  fading  experiment  in  practice  is  so-called  Rayleigh  fading  where  the 
envelope  ‘V2E/T  Is  distributed  In  accordance  with  the  Rayleigh  distribution.  For  that 
case,  equations  (4a)  and  (4b)  become: 


(noncoherent  case) 


and 


^E/No 
2  +  E/No 


(coherent  case) 


respectively,  where  E  is  the  average  value  of  the  random  variable  E. 


(5a) 


(5b) 


It  is  important  to  realize  that  equations  (4)  and  (5)  above  do  not  apply  to  situations 
where  the  envelope  fluctuates  during  one  pulse  period.  Specifically,  the  above  results 
do  not  apply  to  communications  channels  where  the  fading  rate  is  of  the  same  order 
or  faster  than  the  keying  rate. 

Mathematical  expressions  have  been  derived the. solutions  of  which  give  the 
optimum  receiver  for  a  fast  fading  case.  Unfortunately,  these  equations  have  not  been 
solved  in  general  nor  has  the  minimum  probability  of  error  been  estimated.  (Priced) 
has  derived  the  minimum  probability  of  error  for  on-off  keying  with  fast  fading  but 
not  for  FSK. )  The  aim  of  this  report  is  to  evaluate  the  performance  of  a  particular 
(sub-optimum)  receiver  for  this  situation.  The  reason  for  the  choice  of  the  receiver 
chosen  will  be  presented  later. 


2.  MATHEMATICAL  MODEL 


In  the  model  to  be  considered  it  is  assumed  that  the  receiver  has  as  its  input 
waveform  a  narrow  band  Gaussian  signal  centered  either  at  frequency  u>0  or  coi  which 
is  also  corrupted  by  additive  Gaussian  white  noise  (again  with  zero  mean  and  power 
spectral  density  Sq  (co)  =  No/2  for  all  u>).  In  particular,  the  received  waveform  r(t)  is 
given  as 


r(t) 


no  « 
or 


+  n(t> 


0<Lt<  T 


(6) 
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where 


(a)  nQ  (t),  (t)  and  n(t)  are  independent  Gaussian  processes  all  with  zero  mean 

<b)  E  Jn(t)n(t-  r)f  =  ^  fi(T) 

(c)  E  jnQ(t)no  (t  -  r) }  =  Rq  (r) 

E  |n1(t)n1  (t  -  r) }  =  R1  (r) 

(d)  The  a  priori  probabilities  of  nj^  (t)  and  nQ  (t)  are  equal. 

The  receiver  must  operate  on  the  received  waveform  r(t)  during  the  interval 
(0,  T)and  make  a  decision  whether  nQ(t)  or  n1(t)  was  present  during  that  interval. 

The  physical  reasoning  for  such  a  mathematical  model  is  that  when  a  sinusoid  is 
transmitted  over  a  scatter  communications  channel,  it  travels  over  many  different  paths, 
each  path  introducing  amplitude,  phase  and  perhaps  frequency  changes.  The  sum  of  the 
signals  from  all  these  paths  then  (from  Central  Limit  Theorem  arguments)  can  be 
considered  as  a  narrow  band  Gaussian  process  with  center  frequency  given  by  the  fre  ¬ 
quency  of  the  transmitted  carrier.  An  artificial  channel  which  was  constructed  and  ex¬ 
hibits  this  type  of  perturbation  was  the  Needles  belt (6). 

An  intuitive  argument  might  suggest  that  a  natural  method  for  deciding  between  the 
two  noise  sources  is  to  estimate  the  energy  in  the  received  waveform  in  the  two  nar¬ 
row  frequency  bands  centered  at  u>0  and  o)i  and  then  choose  that  frequency  having  the 
largest  energy.  A  breadboard  simulation  of  such  a  system  was  reported  by  Cossette 
and  Wolf  (7).  Theoretical  analyses  of  such  a  system  have  been  made  by  Jacobs  (8)  (for 
one  type  of  spectrum)  and  by  Kobos  and  Meyer  (9). 

Another  sub-optimum  receiver  is  suggested  by  the  following  ideas.  Note  that  the 
correlator  receiver  given  in  Figure  1  cross -correlated  the  received  waveform  with 
stored  replicas  of  the  transmitted  signals.  Thus  it  would  seem  that  a  logical  design 
for  a  receiver  would  be  to  obtain  estimates  of  the  two  noise  waveforms  no(t)  and  nj(t) 
by  passing  the  received  signal  through  filters  centered  at  w0  and  coj  and  then  cross- 
correlating  these  estimates  with  the  received  waveform  r(t).  Such  a  receiver,  called 
a  cross -correlator  receiver,  is  shown  in  Figure  3.  This  receiver  is  analyzed  in  the 
remaining  sections  of  this  report. 

It  should  be  clearly  understood  that  no  claim  for  optimality  is  made  for  this  re¬ 
ceiver.  However,  it  is  interesting  to  note  that  one  form  for  the  block  diagram  of  the 
optimum  receiver  (which  leads  to  minimum  probability  of  error)  is  similar  to  this 
receiver(10).  The  optimum  receiver,  however,  utilizes  time  varying  filters  in  place 
of  narrow  band  filters.  The  time-varying  impulse  responses  for  these  filters  are  not 
known,  but  are  only  known  to  be  the  solutions  to  certain  integral  equations.  Further¬ 
more,  since  time  varying  filters  may  be  difficult  to  build,  there  appears  to  be  ade¬ 
quate  justification  for  considering  this  simpler,  but  sub-optimum,  receiver. 
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Figure  3.  Cross-Correlator  Sub-Optimum  Receiver  for  Noise-Like  FSK  Reception 


3.  MATHEMATICAL  ANALYSIS 

In  the  following  analysis  specific  forms  are  assumed  for  the  transfer  functions  of 
the  filters  in  Figure  3,  as  well  as  for  the  autocorrelation  functions  of  the  noise  pro¬ 
cesses.  Although  the  problem  could  be  analyzed  in  greater  generality,  this  approach 
was  taken  since  a  more  general  analysis  would  be  more  complicated  and  the  salient 
points  are  illustrated  in  the  analysis  which  follows. 

Specifically,  we  assume  that  the  noise  processes  no(t)  and  ni(t)  have  autocor¬ 
relation  functions  identical  to  those  which  would  result  by  passing  white  noise  through 
single-tuned,  high  Q,  RLC  filters  centered  at  0)0  and  Wj  respectively.  Thus,  if  we 
write  no (t).  and  ni(t)  as 


n.(t)  =  Xj(t)  cos  Wj  t  +  y^ft)  sin  0^(1) 

i  =  0. 

1 

(7) 

_ CC  ItI 

Then  R  -R  (r)  =  e  m 

xixi  TO 

1  =  0, 

1 

(8) 

o 

n 

•— » 

tf* 

o 

ii  • 

•H 

1 

0) 
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Note  that  "S”  Is  the  power  In  xi(t)  and  yj(t).  Aa  a  consequence  of  equation  (9),  however, 
it  is  also  the  power  in  ni(t).  Furthermore,  the  impulse  response  of  the  receiver  fil¬ 
ters  are: 


Narrow  band  filter  centered  at 

o 


ho<x)  = 


cos  a>Qx 


x>oi 

x<0) 


(10) 


Narrow  band  filter  centered  at 


e  cos  ojjX  x  >  0 
0  x  <  0 

That  is,  these  receiver  filters  are  just  RLC,  high  Q  filters  centered  at  u>0  and  u>i, 
respectively.  Furthermore,  we  assume  that  the  center  frequencies  o>0  and  coj  are 
separated  far  enough,  so  that  there  is  no  output  of  filter  centered  at  wi  due  to  the  nar¬ 
row  band  noise  n<,(t)  (centered  at  <jj0)  and  there  is  no  output  of  filter  centered  at  cu0 
due  to  the  narrow  band  noise  nj(t)  (centered  at  o>i).  Of  course,  both  filters  have  out¬ 
puts  due  to  the  additive  white  noise  n(t). 

Two  slightly  different  versions  of  the  problem  will  be  considered.  They  are  de¬ 
scribed  below  and  referred  to  as  case  1  and  case  2.  In  both  situations  we  will  assume 
that  the  received  waveform  actually  contained  the  narrow  band  noise  no(t)  plus  the 
additive  Gaussian  white  noise  n(t).  We  then  calculate  the  probability  of  error  as  the 
probability  that  the  receiver  decides  that  ni(t)  was  present.  This,  of  course,  is  just 
one  type  of  error  that  the  receiver  could  make  but  due  to  the  symmetry  of  the  prob¬ 
lem,  this  error  probability  is  equal  to  the  overall  error  probability. 


4.  CASE  1 

Consider  the  receiver  shown  in  Figure  4,  where  the  received  waveform  is  given 
as 


r(t)  =  nQ(t)  +  n(t).  0<t<.T 

The  probability  of  error  for  this  receiver  is  then 


(12)  * 


Pe  =  Pr[Q<°]  'PrlQ^Qj. 


(13) 
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Figure  4.  Correlator  Receiver  for  Case  1. 

As  in  equation  (7),  we  write 

no(t)  =  xo(t)  cos  (i)Qt  +  yQ(t)  sin  a>Qt  (14) 

The  switch  in  the  front  end  of  the  receiver,  which  closes  at  t  =  0,  is  included  to  speci¬ 
fically  indicate  that  the  receiver  filters  are  inert  (have  no  initial  conditions)  at  the 
start  of  the  keying  interval.  (In  practice  the  energy  in  these  filters  would  have  to  be 
dumped  at  the  end  of  every  pulse  period. ) 

To  calculate  the  output  voltage  Q0,  we  write  the  white  noise  n(t)  as* 

n(t)  =  xn(t)  cos  o>0  t  +  yn(t)  sin  u)Qt  (15) 

where  E[xn(t)xn  (t  -  T)  ]  =  E[yn(t)  yQ(t  -  T)]  =  ~~  6  (T)  (16) 

and  E  [xn(t)  yn(t  -  r)]  =o  (17) 


♦The  complications  of  writing  white  noise  in  the  form  given  in  equation  (15)  are  dis¬ 
cussed  by  Helstromf11). 
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Since  r^ft)  and  n(t)  are  statistically  independent  processes,  we  then  have 


r(t)  -  xr(t)  cos  wQt  +  yr  (t)  sin  a>0t 

where 


(18) 


E[  xr(t )  xr  (t  - 


T’ 


]  =  &HTU  >£- 


8  (r>  =  E[yr(t)yr  (t  -  t)]. 


and 


(19) 


E.xr(t)  yr  (t-r)]  =o 


(20) 


Now 


T  (  r  -  Of(t  -  77 ) 

Q  f  <  |  e  cos  o?c(t  -  77 )[  xr(n )  cos  ajQr}  +  yr  (77 )  sin 

0*0 

[xr(t)  cos  o>Qt  +  yr(t)  sin  o)Qt]  dt 


Q>ori]  d77j 


(21) 


But  if  f(x,  y)  =  f(y,  xj,  then 

T 


;r  v  1  /T  T 

/  [  J  f(x,  y)  dx]  dy  =  £  j  /  f(x,  y)  dx  dy 


(22) 


o  o 

so  that  Q  can  be  written  as 
o 


o  o 


T  T  -a/(t  -77) 

%=  I  f  fe  cos  wo(t -77)  [xr(T7)  coswo77  +yr  (77)  sina>or7j 

o  o 


[xr<t)  cosa7Qt  +  yr(t)  sin  0Jot]  dt  <377  (23) 

After  expanding  all  trigonometric  products  in  sums  of  trigonometric  functions,  we  can 
ignore  all  terms  which  involve  cos  w0t,  sin  u>0t,  cos  U)07) ,  or  sinco077  ,  since  their 
contributions  to  Q&  will  be  negligible  after  performing  the  integration.  We  are  then  left 
with  the  expression 

T  T-  alt  -  77 1 

Q0=|  /  /e  1  '[xr(t)xr(77)  +yr(t)yr(77)]  dtd77  (24) 

o  o 
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Let  us  now  consider  the  orthonormal  set  of  functions  0i<t),i  -  1,  2,  .  .  .  which 
are  solutions  to  the  integral  equation 


where 


T  i 

/  e“alt  "  Tl  ^(r)  dr  =  Aj  <^(t)  o<t<.T 


f  0±<t)  0  (t)  dt  =  6 


If  we  expand  xr(t)  and  yr(t)  in  terms  of  these  functions  $j(t)  to  obtain 


(25) 


(26) 


00 


V‘>  = 


no 


yr£)  =  S  yj0j  (t) 


j=i 


(27) 


(23) 


then  it  Is  easy  to  show  that 


E[x.x.]  =  (SX.  +  1*2-  )  5y  =  E  [ypr.]  <29 

We  can  thus  rewrite  Q0  in  terms  of  the  coefficients  x. ,  yj  and  the  functions  (t) 
as  3  3 

%  =  i  /  /  e"  a ^  ^  (xj  +yj)  ^3 ( 77 +yk)  d7?  (3° 


But  making  use  of  Equations  (25)  and  (26),  this  becomes 


OO  00 

1  V  •.  J  .  2  t  * 


Q0"  8  ^ 


Z  S 

j-i  ’  1  1  j=i 


where 


1  /  2  ,  2.x 

cr  s  ^  j+y  j)xj 


(31) 


(32) 
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Since  all  processes  under  consideration  are  Gaussian,  xj  and  yj  are  Gaussian  random 
variables  so  that  e  j  is  chi-square  distributed  with  mean  value 

E[cj]  =  i  Aj(SXj+  If-,.  (33, 

Thus  the  probability  density  function  of  e  j  is  given  as 


i  _€  J*- \ 

P(€j)  =  —  e  J  J  >o. 


and  its  characteristic  function  Me  ^  (jv)  is 


(34) 


(35) 


Since  the  e  j  are  all  statistically  independent,  the  characteristic  function  of  Q0  is  then 


MQ  (jv)  =  it  - - -  (36) 

j  (1  -  jre  ^ 

where  €  j  is  given  in  Equation  (33). 

Let  us  now  concentrate  on  calculating  the  statistics  of  Qj.  Since  the  center  fre¬ 
quency  of  this  filter  is  wi,  it  is  convenient  to  write  the  noise  as 


n(t)  =  x’n(t)  cos  Wjt  +  y*n(t)  sin  Wj  t  (37) 

so  that  the  input  to  the  filter  becomes 

r(t)  =  xQ(t)  cod  w0t  +  y0(t)  sin  wQt  4  x»n(t)  cos.Wjt  +  y'n(t)  sin  (38) 

The  output  Q1  is  then 

Ql=  /  e"a(t  cos  wx(t  -7j)  [r(Tj)]  drjj  r(t)  dt  (39) 
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where  r(t)  is  as  given  in  Equation  38.  Substituting  Equation  (38)  into  (39)  (twice)  and 
ignoring  all  terms  winch  have  sinusoidal  variations,  results  in 


T  t 

J  /  /  <s"a(t " 17  ’  (x'n(7)  >  x'n(t;  +  ya(7, )  yn(t)]  an  dt 


Q,=  5  f  f  [*•  01)**  w+jr*  (i))y*  «]  dn  dt 

Expanding  x'n(t)  and  y’n(t)  in  terms  of  the  functions  0j(t)  previously  defined,  as 


x'n(t*  ’  Z  x'j 
00 

y,nw  =  A  y,j  ^ 


and  substituting  into  Equation  (41)  yields 


<WJ|Xj  l  fry2  ♦  *•/] 


By  a  similar  set  of  steps  to  that  which  led  to  Equation  (36)  we  obtain 


Mq  (jv) 


j  (1-Jve'j) 


where 


f  NoXj 
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We  are  now  In  a  position  to  calculate  the  probability  of  error.  Returning  to  Equa¬ 
tion  (13)  we  see  that  the  probability  of  error  is  given  as 


oo  r  ao 

Pe  =  Pr  [  <3X  >  Qj  -  /  Pfi0>  If  Pl^dQ 

0  LS> 

But  the  probability  density  functions  p(Qj)  and  the  characteristic  functions  Mqj  ()v)  are 
related  by  the  equation 


i  r°  -JvQ.  ,  +)VQ 

p«i>  =  W  L  Mo  S’)®  *>--b  L  m*q,  flv)e  dv 


Substituting  one  form  of  Equation  (48)  into  Equation  (47)  we  obtain 


Pe  =  /  P«„) 


[£[- 


Cv)  e 


+jvQ- 


vjdQx] 


0.0  r  4*o 

Pe  =  /p«o)[2f/„  M\ 


+j  vQ, 


if  Re  Ov)<o. 


Performing  the  integration  with  respect  to  Qq  yields 

+»  mq  Ov)  M*q  flv)  dv 

- - - V-1 -  ifEeflv)<o  (51) 

Finally,  we  can  rewrite  Equation  (51)  in  terms  of  a  complex  variable  s,  as 

M  P)M.  <8, 

to'-m  l  s  -  1)8  <52' 

where  C  is  a  contour  along  the  imaginary  axis  but  to  the  left  of  all  the  poles  on  this 
axis. 
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From  Equations  (36)  and  (45)  we  have  that 


u  /c\  _  lim  N  1 


and 


lim  N 

n— «  rr 


j-i  (1  - s  c ? 


(53) 


(54) 


where 


and 


V  i  xi  <sxi + 


£'i  =  #Vo 


(55) 


(56) 


Thu3  we  can  write  (assuming  that  we  can  take  the  limit  after  integrating) 


Pe  = 


lim 

N — ►°° 


ds 


f _ 

27rj  csfi  (i-s  rp  a  +  s  r3) 


(57) 


i=l 


Defining  the  pole  locations  Sj  and  S'j  as 


S. 

3  r. 


xj  (sxj +  ir  > 


(58) 


S’.  = 


3  €», 


8 


(59) 


we  then  have 


Pe  - 


lim  (-1) 
N — ►*> 


N  -f-  1 


TT  si  s’i  ds 

i=l 


)  f  1=1 _ 

2iri  6  s-fi-  (S-Sj)  (S+S’j) 

i=l 


(60) 
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In  order  to  perform  the  Integration  indicated  in  Equation  (60)  we  note  the  pole  plot 
shown  in  Figure  5.  (The  eigenvalues,  \lt  have  been  ordered  3uch  that  A*  >A->>Ao.  \ 

Let  us  write  Equation  (60)  as  1  ^  J 


Pe  -  2?f  X  fN  * 


Then  from  the  residue  theorem,  the  probability  oi  error  is  the  sum  of  the  residues  of 
fjj(s)  for  the  left  half  plane  poles.  The  residue  of  f^(s)  at  S  =  -S'j  is  given  as: 

TT  i  i 

'  =-S-  =  (-1)N+1  - 55 - - jr -  (62) 

1  <-Sj>  TJ  (-1)  (S:j  +  Sj)  jy  (S'j  -  S'j) 

i=l  i==l 


residue  of  fN(s)  at  S  =  -S'^ 


the  probability  of  error  i3  then 


Finally,  substituting  the  values  for  Sj  and  S’j  given  in  Equations  (58)  and  (59)  yields 
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Figure  5.  Foie  Plot  and  Path  of  Contour  Integration  for  Equation  (60) 

Equation  (65)  is  the  desired  result.  Note  that  this  equation  gives  the  probability  of 
error  in  terms  of  eigenvalues  Aj  of  the  integral  equations  given  by  Equation  (25)  and 
also  the  ratio  S/No  where  S  is  the  power  in  the  "signal",  n^)  and  No  is  the  noise  per 
cycles/second  of  bandwidth  for  the  additive  white  noise.  A  computer  was  used  to  cal¬ 
culate  Pe  for  various  values  of  ( aT)  as  ST/N0  varied  over  a  range  of  values.  These 
results  are  discussed  in  a  later  section. 


5.  CASE  2 

The  second  case  differs  from  the  first  in  the  way  in  which  the  additive  white  noise 
is  treated.  Now  it  is  assumed  that  the  cross -correlations  are  each  preceded  by  ad¬ 
ditional  RLC  filters  which  are  not  dumped  at  the  end  of  each  pulse  period.  Furthermore, 
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the  spectrum  of  the  noise  no(t)  is  assumed  to  be  wide  compared  to  the  bandwidth  of 
the  RLC  filter  centered  at  co0  so  that  it  can  be  considered  as  a  white  noise  input.  How¬ 
ever,  the  separation  between  0>o  and  a is  assumed  wide  enough  so  that  can  be 
ignored  as  an  input  to  the  filters  centered  at  o>]_.  The  actual  receiver  structure  and  the 
circuit  whi  ^h  is  claimed  as  equivalent  for  calculation  purposes  are  shown  in  Figure 
(6a)  and  (6b)  respectively. 

The  analysis  for  Case  2  is  performed  in  a  similar  fashion  to  that  of  Case  1.  The 
resultant  expression  for  probability  of  error  is 


where 


ml  power  +  noise  power 
noise  power 


and  again  Xj  are  the  eigenvalues  of  the  integral  equation  given  in  Equation  (25).  Com¬ 
puter  results  are  also  given  for  Case  2  in  the  next  section. 


6.  COMPUTER  EVALUATION  OF  ERROR  PROBABILITIES 

In  order  to  compare  Cases  1  and  2,  a  common  set  of  parameters  must  be  defined. 
In  the  introduction,  it  was  seen  that  error  probability  expressions  were  always  ex¬ 
pressed  in  terms  of  the  ratio  E/No  where 


E 

No 


energy  in  the  signal  for  one  pulse  period 
noise  power  density 


(68) 


This  ratio  will  also  be  used  for  the  two  cases  considered  here.  The  second  parameter 
was  chosen  as  the  dimensionless  "time-bandwidth"  product,  0,  defined  as 


aT 

2 


(69) 


Furthermore,  it  is  easily  shown  that  the  expressions  for  error  probability  depend  only 
on  the  parameter  /9  and  not  on  the  individual  values  of  a  and  T.  Thus  we  can  arbi¬ 
trarily  set  a  equal  to  10  so  that 

/3  =  5T 

In  terms  of  these  two  parameters  0  and  E/N0)  the  two  expressions  for  probability 
of  error  become 

Case  1 


A  slightly  different  parameter  that  could  be  used  in  describing  the  performance 
of  the  system  is  the  one-sided  three  dfc  bandwidth  "Bn  given  as: 

B  =  =  -pjr  (cycles/ second)  (72) 

In  terms  of  this  parameter,  the  two  probabilities  of  error  become: 

Case  1: 
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where 

g 

^  =  1  +  N  (2BT) 
o 

The  performance  curves  to  be  presented  give  the  probability  of  error  versus  E/N0 
(measured  in  db. )  for  various  values  of  BT. 


7.  COMPUTATION  OF  EIGENVALUES 

The  next  step  is  the  computation  of  the  eigenvalues  (Xj’s)  to  be  substituted  into 
Equations  (70)  and  (71).  These  eigenvalues  which  satisfy  the  integral  equation  given  by 
Equation  (25)  have  been  shown(12T  to  be  related  to  the  non-negative  roots  of  the  trans¬ 
cendental  equations 


tan  z  =  j3/z  (75a) 


cot  z  =  -P/z  (75b) 

As  can  be  seen  from  Figure  7,  the  smallest  root,  Zj,  is  a  solution  of  (75a),  the 
next  smallest  root,  Z2,  is  a  solution  of  (75b),  etc. ,  with  the  roots  alternating  between 
the  two  equations.  The  derived  eigenvalues  are  related  to  the  roots  by  the  equation: 


\  -  h  * 

i 

P  +  z". 


Two  checks  that  were  used  in  the  computation  of  these  eigenvalues  were: 


(a)  The  sum  of  the  eigenvalues  should  equal  T 

th 

(b)  The  n  root  of  Equation  (75a)  should  fall  in  the  interval  (d-1)7T<  z  <. 
(n-l/2)  n  while  the  n1*1  root  of  Equation  (75b)  should  fall  in  the  interval  (n-l/2)?r<z 
^ vcn .Furthermore,  as  n  increases,  the  root  occurs  very  close  to  the  lover  limit  of 
allowed  values. 
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8.  COMPUTATION  OF  ERROR  PROBABILITIES 

The  mathematical  forms  of  the  computation  [Equations  (70),  (71),  (75),  (76)]  were 
programmed  in  FORTRAN  IV.  The  programs  for  cases  1  and  2  are  given  in  the  Ap¬ 
pendix.  BT  products  of  0. 5,  1.0,  2.0,  3.0,  4.0,  and  5,  0  were  evaluated  for  E/N0 
ratios  of  6  db  through  24  db  in  increments  of  two  db.  Tabulations  of  these  outcomes 
follow  with  their  corresponding  graphs.  (Figures  S  through  11.) 


9.  CALCULATIONS  FOR  OTHER  FILTERS 

The  techniques  outlined  in  this  paper  can  be  used  for  filter  forms  other  than  the 
XILC.  in  particular,  if  the  noise  processes  n0(t)  and  ni(t),  as  given  in  equation  (7),  are 
such  that 


Sx.x.  <"> 
1  1 


s 


Vi 


^o,  |o)|>  27TB 

( "2B"  »  |t*>|<_27TB 
I  O,  |  Ct)  j  >  27TB 


(77) 


(78) 


and  if  the  receiver  filters  given  in  Figures  4  and  6  have  transfer  functions  which  have 
unity  gain  over  a  pass-band  of  bandwidth  2B  (cycles/second)  centered  at  u)p  and  o)i, 
and  zero  gain  elsewhere,  then  the  error  probabilities  given  in  Equations  (73)  and  (74) 
still  apply  if  the  Aj  are  the  solutions  to  the  integral  equation 


fT  si2f&rJ  <T) dT  =  \  V>  <79> 

o 

Eigenvalues  for  this  integral  equation  with  a  different  normalization  are  given  by 
Jacobs 


To  arrive  at  an  appropriate  normalization  consider  the  following.  In  the  cases 
first  examined  a  one-sided  filter  bandwidth  was  used  exclusively,  that  is, 

3  r:  _2L 

27T 

For  the  case  of  the  flat  spectrum,  the  double-sided  bandwidth  is  used 


B’ 


QL 

V 


(80) 


Figure  8 


so  that  when  we  set  Of  =  10,  then  B'  =  10/tt.  The  flat  spectrum  eigenvalues,  which  art 
called  Xj's,  are  such  that 


(81) 
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PROBABILITY  OF  ERROR  TABULATION 
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PROBABILITY  OF  ERROR  TABULATION 


Figure  12 


CASE  2:  FLAT  SPECTRUM 
EIGENVALUES 


E/ N.  (db) 


PROBABILITY  OP  EP.ROR  TABULATION 
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Thus,  to  get  correct  results,  we  nust  divide  the  flat  spectrum  eigenvalues  by  B\  or 


i 


\=  <\V  To  <83> 

Note  that  when  the  flat  spectrum  eigenvalues  are  used,  the  normalization  in  Equation 
(83)  will  not  make  any  difference  3ince  the  eigenvalues  always  occur  as  X2j/X2j,  there¬ 
by  canceling  the  normalization.  It  will,  however,  make  a  difference  in  case  one. 


These  eigenvalues  normalized  such  that  )  A  -  T  are  listed  in  Table  ?  (Appendix). 

T*  i 

The  results  obtained  using  these  eigenvalues  for  the  two  cases  are  given  in  Figures  13 
and  15. 


10.  CONCLUSIONS 

This  report  considers  a  sub-optimum  receiver  for  binary  frequency  shift  keyed 
(FSK)  transmission  which  experiences  fast  fading  and  additive  Gaussian  white  noise. 

The  received  signal  is  modeled  as  one  of  two  narrow  band  Gaussian  processes  (whose 
center  frequencies  correspond  to  the  frequencies  used  in  the  binary  FSK  transmission) 
corrupted  by  additive  Gaussian  white  noise.  The  half -bandwidth  of  each  narrow  band 
process  is  B  (ops)  and  the  pulse  duration  is  T  (seconds). 

A  receiver  is  considered  which  passes  the  received  waveforms  through  one  of  two 
narrow  band  filters  (with  center  frequencies  as  above)  aDd  cross -correlates  the  outputs 
of  these  filters  with  the  received  waveform.  The  receiver  makes  a  decision  every  T 
seconds  as  to  which  frequency  was  transmitted  by  choosing  that  frequency  correspond¬ 
ing  to  the  larger  value  of  the  crosa -correlation.  Two  slightly  different  forms  of  this 
receiver  are  analyzed  and  curves  for  the  probability  of  error  versus  the  ratio  of  energy 
per  bit  to  noise  power  density  (E/N0)  for  each  receiver  are  presented  with  the  product 
BT  as  a  parameter.  Two  different  spectra  for  the  narrow  band  processes  are  studied. 

The  following  conclusions  were  derived  from  the  results.  First,  increasing  BT 
does  not  give  a  constant  improvement  in  probability  of  error  for  a  fixed  value  of  E/NQ. 
In  fact,  the  rate  of  improvement  diminishes  as  BT  increases.  It  can  further  be  noted 
that  for  E/Nq  greater  than  about  13  db  a  higher  BT  corresponds  to  a  better  system  with 
the  opposite  effect  for  lower  E/N0.  The  results  for  BT  -  0. 5,  1. ,  and  2.  (for  a  parti¬ 
cular  spectrum)  appear  to  be  in  essential  agreement  with  those  obtained  experimentally 
by  Cossette  and  Wolf^)  who  used  a  slightly  different  receiver. 

The  results  of  this  report  would  have  practical  importance  if  it  were  found  that  the 
receivers  analyzed  were  easier  to  construct  than  the  conventional  receiver  for  this  type 
of  modulator  and  channel  behavior.  However,  the  practical  aspects  of  the  receiver  de¬ 
sign  lay  outside  the  realm  of  this  study. 
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DETECTION  EXTENSION 

CORRELATION  RECEIVER 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


C 

C 

C 

c 

c 

c 

c 

c 

c 


13 

14 

15 

16 

17 

18 

19 

20  401 


21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32  633 


ENERGY  DETECTION  EXTENSION 
CROSS  CORRELATION  RECEIVER 
(CASE  1) 

RLC  EIGENVALUES 

PE  MAGNITUDE,  GT,  SMALLEST  EIGEN.  N  =  30 
SINGLE  PRECISION 

PE  VS  EN0DB  =  6,  8,  10,  12,  14,  16,  18,  20,  22,  24 
FIND  EIGENVALUES 
REAL  LAMBDA 

DIMENSION  LAMBDA  (65),  R(16),  AT(20 

DATA  (R(K),  K  =  1,  10)/.  6,-8,  1.,  1.2,  1.4,  1.6,  1.8, 

2.,  2.2,  2.4/ 

DATA  (AT (KK),  KK  =  1,  6)/.  5,  1.,  2.,  3.,  4.,  5./ 

PI  =  3. 14159265 
DO  993  KK  =  1,  6 
BT  =  AT  (KK) 

BETA  =  PI  *  BT 
T  =  0. 62831853*BT 
PRINT  401,  T,  BETA,  BT 

FORMAT  <2X,  3H  T  =  E14.  8,  2X,  6H  BETA  =  E14.  8, 

2X,  4H  BT  =  E14.8) 

N  =  30 

ACCURA  =  1.  E  -  7 
DELTA  =  0.5 
Z  =  0. 

Cl  =  T  *  BETA 
C2  =  BETA  *  BETA 
K2SIGN  =  1 
K1SIGN  =  K2SIGN 
I  =  0 

NNMAX  =  N/2 
PRINT  633 

FORMAT  (3X,  2H,  I,  14X,  2H  Z,  14X,  10K  LAMBDA (I)) 


33 

34 

35 

36 

37 

38 

39 

40 

41 

42 


3 
1 
2 

4 

6 

8 

9 


DO  30  NN  =  1,  NNMAX 
DEL  =  DELTA 
Z  =  Z  +  DEL 

FI  =  BETA/Z  -  SIN(Z)/ COS(Z) 

IF (Fl)4,  4,  6 
LlSIGN  =  0 
GO  TO  8 
LlSIGN  =  1 

IF  (K1SIGN  -  LlSIGN)  9,  1,  9 
IF(ABS(DEL)  -  ACCURA)11,  11,  10 
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DETECTION  EXTENSION 


CORRELATION  RECEIVER 


43 

10 

K1SIGN  =  L1SIGN 

44 

DEL  =  -0. 1*DEL 

45 

GO  TO  2 

46 

11 

1=1  +  1 

47 

LAMBDA  (I)  =  C1/(C2  +  Z*Z) 

48 

152 

CONTINUE 

49 

645 

PRINT  5050,  I,  Z,  LAMBDA (I) 

50 

5050 

FORMAT  (2X,  15,  2E22.8) 

51 

644 

CONTINUE 

52 

K1SIGN  =  L1SIGN 

53 

23 

DEL  =  DELTA 

54 

21 

Z  =  Z+DEL 

55 

22 

F2  =  (COS(Z)/SIN(Z)  )  +  (BETA/Z) 

56 

IF  (F2)  24,  24,  26 

57 

24 

L2SJGN  =  0 

58 

GO  TO  28 

59 

26 

L2SIGN  =  1 

60 

28 

IF (K2 SIGN  -  L2SIGN)29,  21,  29 

61 

29 

IF(ABS(DEL)  -  ACCURA)41,  41,  20 

62 

20 

K2SIGN  =  L2SIGN 

63 

DEL  =  -0. 1*DEL 

64 

GO  TO  22 

65 

41 

1  =  1  +  1 

66 

LAMBDA (I)  =  C1/{C2  +  Z  *  Z) 

67 

1645 

PRINT  5050,  I,  Z,  LAMBDA (I) 

68 

1644 

CONTINUE 

69 

30 

K2SIGN  =  L2SIGN 

70 

SSUM  =  0,  0 

71 

DO  4049  L  =  1,  N 

72 

4049 

SSUM  =  SSUM  +  LAMBDA  (L) 

73 

PRINT  9099,  SSUM 

74 

9099 

FORMAT  (2X,  18H  LAMBDA  SUM  EQUALS,  E20.  8) 

75 

D  =  N  -  1 

76 

TEST  =  D*PI/2 

77 

PRINT  2020,  TEST 

78 

2020 

FORMAT  (2X,  6H  TEST  =  E20.  8) 

79 

C 

CALCULATE  PE 

80 

PRINT  4036 

81 

4036 

FORMAT  (2X,  2H  I,  14X,  5H  GLOB) 

82 

DO  501  K  =  1,  10 

83 

EN0DB  =  R(K)*1C 

84 

EN0  =  10.  **R(K) 

85 

SUM  =  0. 

86 

DO  900  I  =  1,  N 

87 

GLOB  =  1.  E25 

88 

ACON  =  1.  0/ LAMBDA  (I) 

35 
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89 

CHI  =  LAMBDA  (I)/BETA 

90 

SI  =  EN0*CHI 

91 

W  =  1.  +5.  *SI 

92 

WT  =  l./W 

93 

PROD  =  1. 

94 

DO  200  J  =  1,  N 

95 

IF (J -I)  69,  200,  69 

96 

69 

X  =  LAMBDA  (J)*ACON 

97 

CHJ  =  LAMBDA(J)/BETA 

98 

SJ  =  EN0*CHJ 

99 

V  =  1.  +10.  *SJ 

100 

VT  =  V*X 

101 

VTH  =  1.  +VT 

102 

XX  =  1.  -X 

103 

PROD  =  PROD*VTH*XX 

104 

200 

CONTINUE 

105 

FKTR  =  l./PROD 

106 

GLOB  =  GLOB*WT*FKTR 

107 

999 

PRINT  9090,  I,  GLOB 

108 

9090 

FORMAT  <2X,  15,  E22.  8) 

109 

900 

SUM  =  SUM  +  GLOB 

110 

899 

PE  =  0,  5*SUM*1.  E-25 

111 

100 

PRINT  111,  PE,  EN0DB,  BETA 

112 

111 

FORMAT  <2X,  4H  PE  =,  E29.  8,  5X,  7H  EN0DB 
E20.  8,  5X,  6H  BETA  =  E20.  8) 

113 

501 

CONTINUE 

114 

993 

CONTINUE 

115 

STOP 

116 

END 
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1 

9 

3 

4 

5 

6 

7 

8 
9 

10 

11 


12 

13 

14 
J.5 
36 

17 

18 

19 

20 
21 


C 

C 

c 

c 

c 

c 

c 

c 


401 


22 

23 

24 

25 
2S 

27 

28 

29 

30 

31 

32 

33  633 


34 

35 

36 

37 

38 

39 

40 

41 

42 

43 


3 
1 
2 

4 

6 

8 

9 


CORRELATION  RECEIVER 

ENERGY  DETECTION  EXTENSION 
CROSS  CORRELATION  RECEIVER 
RLC  EIGENVALUES  (CASE  2) 

PE  MAGNITUDE,  GT,  SMALLEST  EIGEN.  N  »  30 
SINGLE  PRECISION 

PE  VS  EN0DB  =  6,  8,  10,  12,  14,  16,  18,  20v  22,  24 
BT  =  .5,  1.,  2.,  3.,  4. ,  5. 

FIND  EIGENVALUES 
REAL  LAMBDA 

DIMENSION  LAMBDA  (40),  BAMBDA  (40),  R(16),  AT(10) 
DATA(R(K),  K=  1,  10)/.  6,  .8,  1.,  1.2,  1.4,  1.6,  1.8, 
2. ,  2.2,  2.4/ 

DATA  (AT  (KK),  KK  =  1,  6)/.  5,  1.,  2.,  3.,  4.,  5./ 

PI  =  3.14159265 
DO  993  KK  =  1,  6 
BT  =  AT(KK) 

BTT  =  2.  *BT 

BBTT  =  1.  /BTT 

T  =  0.  62831853*BT 

BETA  =  5.  *T 

PRINT  401,  T,  BETA,  BT 

FORMAT  (2X,  3H  T  =  E14.  8,  2X,  6H  BETA  =  E14.  8, 
2X,  4H  BT  =  E14.  8) 

N  =  30 

ACCURA  =  1.  E-7 
DELTA  =  0.  5 
Z  =  0 

Cl  =  T*BETA 
C2  =  BETA*BETA 
K2SIGN  =  1 
K1SIGN  =  K2SIGN 
I  =  0 

NNMAX  =  N/2 
PRINT  633 

FORMAT  (3X,  2H  I,  14X,  2H  Z,  14X,  10H  LAMBDA (I), 
14X,  10H  BAMBDA(I)) 

DO  30  NN  =  1,  NNMAX 
DEL  =  DELTA 
Z  =  Z+DEL 

FI  =  BETA/Z-SIN'Z^OO::.~4 
IF(F1)^ ,  4,  6 
T  i  =  0 
GO  TO  8 
L1SIGN  =  1 

IF  (K1SIGN-  L1SIGN)9,  1,  9 
IF(ABS(DEL)  -  ACCURA)11,  11,  10 


37 
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44 

10 

K1SIGN  =  L1SIGN 

45 

DEL  =  -Q.  1*DEL 

46 

GO  TO  2 

47 

11 

1  =  1  +  1 

48 

LAMBDA  (I)  =  C1/(C2  +  Z*Z) 

49 

BAMBDA(I)  =  LAMBDA(I)*LAMBDA(I) 

50 

152 

CONTINUE 

51 

645 

PRINT  5050,  I,  Z,  LAMBDA (I),  BAMBDA(I) 

52 

5050 

FORMAT  (2X,  15,  3E22.  8) 

53 

644 

CONTINUE 

54 

K1SIGN  =  L1SIGN 

55 

23 

DEL  =  DELTA 

56 

21 

Z  =  Z+DEL 

57 

ejn 

F2  =  (COS(Z)/SIN(Z))  +(BETA/Z) 

58 

IF(F2)24,  24,  26 

59 

24 

L2SIGN  =  0 

60 

GO  TO  28 

61 

26 

L2SIGN  =  3, 

62 

28 

IFfK2SIGN-L2SIGN)29,  21,  29 

63 

29 

IF(ABS(DEL)  -  ACCURA)41,  41,  20 

44 

20 

K2SIGN  =  L2SIGN 

65 

DEL  =  -0.  1*DEL 

66 

GO  TO  22 

67 

41 

1  =  1+1 

68 

LAMBDA  (I)  =  C1/(C2  +  Z*Z) 

69 

BAMBDA(I)  =  LAMBDA  (I)*  LAMBDA  (I) 

70 

1645 

PRINT  5050,  I,  Z,  LAMBDA  (I),  BAMBBA(I) 

71 

1644 

CONTINUE 

72 

30 

K2SIGN  =  L2SIGN 

73 

SSUM  =  0,  0 

74 

DO  4049  L  =  1,  N 

75 

4049 

SSUM  =  SSUM  +  LAMBDA  (L) 

76 

PRINT  9099,  SSUM 

77 

9099 

FORMAT  (2X,  18H  LAMBDA  SUM  EQUAJJS,  E20.  8) 

78 

D  =  N  -  1 

79 

TEST  =  D*PI/2, 

80 

PRINT  2020,  TEST 

81 

2020 

FORMAT  (2X,  6H  TEST  =  *20,  8) 

82 

C 

CALCULATE  PE 

83 

PRINT  4036 

84 

4036 

FORMAT  (2X,  2H  I,  14X,  5H  PROD) 

85 

DO  501  K  =  1,  10 

86 

FNDB  =  R<K)*10. 

G7 

EN  =  10.  **R(K) 

88 

SN  =  BBTT*EN 

89 

SNDB  =  10.  *AL0G10(SN) 

38 
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90 

U  =  1.  +  SN 

91 

UU  =  U-l. 

92 

SUM  =  0. 

93 

DO  900  1  =  1,  N 

34 

999 

PRINT  9090,  I,  PROD 

95 

9090 

FORMAT  (2X,  15,  E22.  8) 

96 

ACON  =  1.  0/ LA  MB  DA  (I) 

97 

PROD  =  1.  0E35 

98 

DO  200  J  =  1,  N 

99 

IF(J-I)  69,  200,  69 

100 

69 

X  =  LAMBDA  <J)*ACON 

101 

Y  =  X*X 

102 

PROD  =  PROD*l.  0/(1.  0  +  Y*(UU-U*Y)> 

103 

200 

CONTINUE 

104 

900 

SUM  =  SUM  +  PROD 

105 

899 

PE  =  1.  0/(1,  0  +  U)*SUM*1.  OE-35 

106 

100 

PRINT  111,  PE,  ENDB,  SNDB 

107 

111 

FORMAT  (2X,  4H  PE  =,  E20.  8,  5X,  6H  ENDB  =,  E20.8 
5X,  6H  SNDB  =,  E20.8) 

108 

501 

CONTINUE 

109 

993 

CONTINUE 

110 

STOP 

111 

END 

39 
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I 

C 

FLAT  SPECTRUM  EIGENVALUES  CASE  1 

9 

REAL  LAMBDA 

3 

DIMENSION  LAMBDA (50,  R(16),  AT(10) 

4 

DATA(R(K),  K  =  1,  10)/.  6,  .8,  1.,  1.2,  1.4,  1.6,  1 
2.,  2.2,  2,  4/ 

5 

PI  =  3.  14159265 

6 

DO  993  KK  =  1,  4 

7 

BT  =  KK 

8 

BETA  =  PI*BT 

9 

BTT  =  2,  ’BT 

10 

BBTT  =  1.  /BTT 

11 

READ  931,  N 

12 

931 

FORMAT  (12) 

13 

READ  246,  (LAMBDA (I),  1=1.  N) 

34 

FORMAT  K 

PRINT  247,  (LAMBDA  (I),  1  =  1,  N* 

10 

247 

format  (E20.  8) 

17 

SSUM  =  0.  0 

18 

DO  4049  L  =  1,  N 

19 

4049 

SSUM  =  SSUM  +  LAMBDA  (L) 

20 

PRINT  9099,  SSUM 

21 

9099 

FORMAT  (2X,  18 H  LAMBDA  SUM  EQUALS,  E20.  8) 

.22 

C 

CALCULATE  PE 

23 

PRINT  4036 

24 

4036 

FORMAT  (2X,  2H  I,  14X,  5H  GLOB) 

25 

DO  501  K  =  1,  10 

26 

EN0DB  =  R(K)*10. 

27 

EN0  =  10  **R(K) 

28 

SUM  =  0 

29 

DO  900  I  =  1,  N 

30 

GLOB  =  1.  E25 

31 

ACON  =  1. 0/ LAMBDA  (I) 

32 

CHI  =  LAMBDA  (I)/BET  A 

33 

SI  =  EN0*CHI 

34 

W  =  1.  +5,  *SI 

35 

WT  -  l./W 

36 

PROD  =  1 

37 

200  J  =  1,  N 

38 

JF(J-1)  63,  200,  69 

3G 

69 

X  =  LAMBDA  (J)*ACON 

40 

CHJ  =  LAMBDA  (J)/DETA 

41 

SJ  =  EN0*CHJ 

42 

V=  1.  +10.  *SJ 

43 

VT  =  V*X 

44 

VTH  -  1.  +VT 

45 

XX  =  1.  -X 

46 

PROD  =  PROD*VTH*XX 

40 
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47 

200 

CONTINUE 

48 

FKTR  =  1,/PROD 

49 

GLOB  =  GLOB  *WT*  FKTR 

50 

999 

PRINT  9090,  I,  GLOB 

51 

9090 

FORMAT  (2X,  15,  E22.  8) 

52 

900 

SUM  =  SUM  +  GLOB 

53 

899 

PE  =  0. 5*SUM*1.  E-25 

54 

100 

PRINT  111,  PE,  EN0DB,  BETA 

55 

111 

FORMAT  (2X,  4H  PE  =,  E20.  8,  5X,  7H  EN0DB  =,  E20  .8 
5X,  6H  BETA  =,  E20.  8) 

5  C 

T  *  \  *5 

CONTINUE 

57 

993 

CONTINUE 

58 

STOP 

59 

END 

41 


1 

t  K0021  1 

04-24-67 

DETECTION  EXTENSION 

CORRELATION  RECEIVER 

'  1 

C 

FLAT  SPECTRUM  EIGENVALUES  CASE  E 

2 

REAL  LAMBDA 

3 

DIMENSION  LAMBDA (50),  R(16),  AT(10) 

•  4 

DATA(R(K),  K=  1,  10)/.  6,  .8,  1.,  1.2,  1.4,  1.6,  1.8, 

■ 

2.  t  2.2,  2.4/ 

5 

PI  **  3. 14  159265 

6 

DO  993  KK=  1,  4 

7 

BT  =  KK 

g 

BTT  -  ■  2.  +J2T 

\  9 

BBTT  =  1,/BTT 

i  10 

READ  931,  N 

11 

931 

FORMAT  (12) 

12 

READ  246,  (LAMBDA  (I),  1  =  1,  N) 

13 

246 

FORMAT  (E10,  8) 

14 

PRINT  247,  (LAMBDA (I),  I  =  1,  N) 

15 

247 

FORMAT  (E20,  8) 

16 

SSUM  =  0.  0 

1  17 

DO  4049  L  -  1,  N 

18 

4049 

SSUM  =  SSUM  +  LAMBDA  (L) 

19 

PRINT  9099,  SSUM 

,  20 

9099 

FORMAT  (2X,  18H  LAMBDA  SUM  EQUALS,  E2Q  8) 

21 

C 

CALCULATE  PE 

22 

PRINT  4036 

1  23 

4036 

FORMAT  (2X,  2H  I,  14X,  5H  PROD) 

1  24 

DO  501  K  =  1,  10 

25 

ENDB  =  R(K)*10 

i  26 

EN  =  1Q.**R(K) 

t  27 

SN  =  BBTT*EN 

28 

SNDB  =  10.  *ALOG10(SN) 

29 

U  =  1.  +  SN 

30 

UU  =  U  -  1 

31 

SUM  =  0. 

32 

DO  900  1=  1,  N 

33 

999 

PRINT  9090,  I,  PROD 

34 

9090 

FORMAT  (2X,  15,  E22.  8) 

35 

ACON  =  1.  0/ LAMBDA (I) 

36 

PROD  =  1.  0E35 

37 

DO  200  J  =  1,  N 

38 

IF(J-I)  69,  200,  69 

39 

69 

X  =  LAMBDA  (J)*ACON 

40 

Y  =  X*X 

.  41 

PROD  =  PROD*l.  0/(1.  0  +  Y*  (UU-U*Y) ) 

42 

200 

CONTINUE 

43 

900 

SUM  =  SUM  +  PROD 

44 

899 

PE  =■  1.  0/(1.  0  +  U)*STJM*1.  0E-35 

45 

i 

t 

100 

PRINT  111,  PE,  ENDB,  SNDB 

42 


O  rfk  *4.  r*. 
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111  FORMAT  (2X,  4H  PE  =,  E20.  8,  5X,  6H  ENDB  =,  E20  .8 

5X,  6H  SNDB  E20.  8) 

501  CONTINUE 

993  CONTINUE 

STOP 
END 


TABLE  #1 


RLC 


LAMBD 
TEST  = 


1415927E  00  BETA  =  0. 15707963E  01  BT  =  0.  50000000E  00 


I 

Z 

LAMBDA  (I) 

1 

0. 10026742E 

01 

0. 14210043E  00 

2 

0. 21924764E 

01 

0.  67838358E-01 

3 

0. 35574032E 

01 

0.  32632144E-01 

4 

0. 50158793E 

01 

0.  .17862602E-01 

5 

0. 65196132E 

01 

0. 10972870E-01 

6 

0. 80467656E 

01 

0.  73415069E-02 

7 

0.  95871783E 

01 

0. 52285754E-02 

8 

0. 11135709E 

02 

0.39019110S-02 

9 

0.  12689531E 

02 

0. 30183328E-02 

10 

0. 14246978E 

02 

0.  24020207E-02 

11 

0.  15807011E 

02 

0. 19557015E-02 

1? 

0. 17368951E 

02 

0.  16225010E-02 

13 

0. 18932335E 

02 

0. 13673573E-02 

14 

0.204  96839E 

02 

0. 11677579E-02 

15 

0.22062227E 

02 

0. 1U087303E-02 

16 

0. 23628326E 

02 

0. 88001218E-03 

17 

0. 25I95006E 

02 

0.  77438333E-03 

18 

0. 26762165E 

02 

0.  68664722E-03 

19 

0. 28329724E 

02 

0.  61298793E-03 

20 

0. 29897621E 

02 

0.  55055324E-03 

21 

0. 31465806E 

02 

0.497 17705E-03 

22 

0. 33034237E 

02 

0. 45119089E-03 

23 

0. 34602883E 

02 

0. 4 112  92 84E -03 

24 

0.36171714E 

02 

0. 37645523E-03 

25 

0.377407082 

02 

0.  34585812E-03 

26 

0. 39309846E 

02 

0.  31884099E-03 

27 

0. 4087911  IE 

02 

0.  29486695E-03 

23 

0.42448488E 

02 

0. 27349601E-03 

29 

0. 44017967E 

02 

C.25436484E-03 

30, 

A  SUM  EQUALS 

0. 45587536E 

02 

0. 30738609E  00 

0.23717133E-03 

0. 45553093E  02 


44 


TABLE  #2 


RLC 


T  =  0.62831853E  00 

BETA  =  0. 31415927E  01 

BT  =  0. 10000000E  01 

I 

Z 

LAMBDA  (I) 

1 

0. 12046393E 

01 

0. 17436293E  00 

2 

0. 24744355E 

01 

0. 12342841E  00 

3 

0. 38287281E 

01 

0.  80473722E-01 

4 

0.52515039E 

01 

0, 52711125E-01 

5 

0. 67204689E 

01 

0. 35867098E-01 

6 

0. 82190773E 

01 

0. 25495331E-01 

7 

0. 97368817E 

01 

0. 18857350E-01 

8 

0. 11267488E 

02 

0. 144265 16E -01 

9 

0. 12806925E 

02 

0. 11351750E-01 

10 

0. 14352654E 

02 

0.  91440984E-02 

11 

0. 15902999E 

02 

0. 75118293E-02 

12 

0. 17456817E 

02 

0.  62741821E-02 

13 

0. 19013308E 

02 

0.  53151622E-02 

14 

0. 20571894E 

02 

0. 45579462E-02 

15 

0. 22132154E 

02 

0. 39501974E-02 

16 

0. 23693767E 

02 

0. 34553553E-02 

17 

0.25256493E 

02 

0. 30473026E-02 

18 

0.26820142E 

02 

0. 27070064E-02 

19 

0. 28384565E 

02 

0. 24203462E-02 

26 

0.29949644E 

02 

0. 21766766E-02 

21 

0.31515233E 

02 

C.  19678546E-02 

22 

0.33081404E 

02 

0. 17875688E-02 

23 

0. 34647943E 

02 

0. 16308684E-02 

24 

0.S6214848E 

02 

0. 14938274E-02 

25 

0.37782071E 

02 

0. 13733013E-02 

26 

0. 39349577E 

02 

0. 12667478E-02 

27 

0.40917333E 

02 

0. 11720940E-02 

28 

0.42485312E 

02 

0. 10876368E-02 

29 

0.44053490E 

02 

0. 10119663E-02 

3C 

0. 45621846E 

02 

0.  94390765E-03 

LAMBDA  SUM  EQUALS 

Q.  60127003E 

00 

TEST  =  0.45553093E  02 


45 


TABLE  #3 


KLC 


T  =  C. 12 56637 IE  01 

BETA  =  0.62831853E  01 

BT  =  0.  20000000E  01 

I 

Z 

LAMBDA  (I) 

1 

0. 13579462E 

01 

0. 19107498E  00 

o 

C4 

0. 27315071E 

01 

0.16820962E  00 

3 

0.41308029E 

01 

0. 13964293E  00 

4 

0. 555S8752E 

01 

0. 11218724E  00 

5 

0. 70136998E 

01 

0. 89045309E-01 

6 

0. 84910380E 

01 

0. 70764979E-01 

7 

0. 99863744E 

01 

0. 56719382E-01 

8 

0. 11495776E 

02 

0. 46003753E-01 

9 

0. 13016102E 

02 

0. 37796957E-01 

10 

0. 14544938E 

02 

0.31452689E-01 

11 

0. 16080456E 

02 

0.  26490293E-01 

12 

0. 17621274E 

02 

0.  22559888E-01 

13 

0. 19166340E 

02 

0. 19407963E-01 

14 

0.20714S50E 

02 

0. 16850112E-01 

15 

0.22266182E 

02 

0. 14751083E-01 

16 

0.23819850E 

02 

0. 1301062  3E-01 

17 

0. 25375488E 

02 

0. 11553654E-01 

13 

0.26932730E 

02 

0. 10323174E-01 

19 

0.  28491389E 

02 

G.  027P5536E-02 

20 

0.30051243E 

02 

0. 83768883E-02 

21 

0.31612128E 

02 

0. 76007364E-02 

22 

0. 33173907E 

02 

0.  69261201E-02 

23 

0.34736465E 

02 

0.  63363137E-02 

24 

0.362 99709E 

02 

0. 58178530E-02 

25 

0. 37863556E 

02 

0. 53598108E-02 

26 

0.39427938E 

02 

0. 49532508E-02 

27 

0.40992796E 

02 

0. 45908133E-02 

28 

0.42558080E 

02 

0. 42663952E-02 

29 

0. 44123745E 

02 

0. 39749042E-02 

30 

0. 456897 54E 

02 

0. 37120660E-02 

LAMBDA  SUM  EQUALS 

0. 11490352E 

01 

TEST  =  0.45553093E  02 


46 


TABLE  #4 


RLC 


T  =  0.  13849556E  01 

BETA  ==  0.  94247780E  01 

BT  =  0.30000000E  0 

I 

Z 

LAMBDA  (I) 

1 

0. 14211365E  01 

0. 19555373E  00 

2 

0. 2848 123 IE  01 

0. 18326397F  00 

3 

0.42856183E  01 

0. 16573185E  00 

4 

0. 57364269E  01 

0. 14593646E  00 

5 

0. 72015146E  01 

0. 12627414E  00 

6 

0. 86804669E  01 

0. 10820817E  00 

7 

0. 10172062E  02 

0. 92384499E-01 

8 

0. 11674738E  02 

0.  78912660E-01 

9 

0. 13186899E  02 

0.  67620420E-01 

10 

0. 1470707 IE  02 

0.  58223103E-01 

11 

0.16233965E  02 

0. 50416849E-01 

12 

0. 177664 93E  02 

0.4392 1874E-01 

13 

0. 19303748E  02 

0. 38497909E-01 

14 

0.2  084498  IE  02 

0. 33945997E-Q1 

15 

0.22389579E  02 

0. 30104549E-01 

16 

0. 23  937036E  02 

0. 26843562E-01 

17 

0.25486935E  02 

0. 24058830E-01 

18 

0.27038932E  02 

0. 21666828E-01 

19 

0.28592740E  02 

0. 19600454E-01 

20 

0.30148120E  02 

0n  17805607E-01 

21 

0.31704873E  02 

0. 16238460E-01 

22 

0, 33262828E  02 

0. 14863337F-01 

23 

0.34821843E  02 

0.13651039E-C1 

24 

0. 36381795E  02 

0. 12577538E-01 

25 

0.37942580E  02 

0. 11622954E-01 

26 

0.39504107E  02 

0. 10770752E-01 

27 

0.41066299E  02 

0. 10007106E-01 

28 

0.42629088E  02 

0.  93203967E-02 

29 

0.44192416E  02 

0. 87008141E-02 

30 

0-.45755230E  02 

0. 81400325E-02 

LAMBDA  SUM  EQUALS 

0. 16448637E 

01 

TEST  =  0.45553093E  02 


47 


TABLE  #5 


BLC 


T  =  0„25132741E  01 

BETA  =  0. 12566371E  02 

BT  =  0. 40000000E 

I 

Z 

LAMBDA  (I) 

1 

0. 14554862E  01 

0.19735248E  00 

2 

0. 29137500E  01 

0.16979597E  00 

3 

0. 43772054E  01 

0. 17835936E  00 

4 

0. 58476472E  01 

0. 16440038E  00 

5 

0.73261580E  01 

0. 14926647E  00 

6 

0. 8813X605E  01 

0. 13406062E  00 

7 

0. 10308560E  02 

0. 11955002E  00 

8 

0. 11811911E  02 

0.10618369E  00 

9 

0.13322568E  02 

0. 94163 125E-01 

10 

0. 1483 980 IE  02 

0. 83522736E-01 

11 

0. 16362873E  02 

0. 74197565E-0I 

12 

0. 17891Q82E  02 

0. 66072043E-01 

13 

(M9423789E  02 

0. 59011389E-01 

14 

0*  20960425E  02 

0. 52879985E-01 

15 

0. 22500489E  02 

0. 47551061E-01 

16 

0. 24043548E  02 

0. 42910981E-01 

17 

0.25589228E  02 

0.38860*28E-01 

18 

0.27137205E  02 

0. 35313949E-01 

19 

0. 28687204E  02 

0.321987.16E-01 

20 

0. 30238985E  02 

0. 29452991E-01 

21 

0.31792343E  02 

0. 27024584E-01 

22 

0, 33347102E  02 

0. 24869413E-01 

23 

0. 34903 106E  02 

0. 22950225E-01 

24 

0. 36460225E  02 

0. 212354  93E-01 

25 

0.38018341E  02 

0. 19698483E-01 

26 

0. 39577354E  02 

0. 1831C469E-01 

27 

0, 41137176E  02 

0. 17070079E-01 

28 

0.42697729E  02 

0. 15942745E-01 

29 

0. 44258 945E  02 

0. 14920263E  01 

30 

0. 458207 63E  02 

0. 13990410E-01 

LAMBDA  SUM  EQUALS 

0. 20911219E  01 

TEST  =  0.45553093E  02 

V 
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TABLE  #6 


RLC 


T  =  0. 31415927E  01 

BETA  =  0. 15707963E  02 

BT  =  0. 50000000E  01 

I 

Z 

LAMBDA  (I) 

i 

0. 14770406E 

01 

0. 19824712E  00 

2 

0.29556073E 

01 

0. 19316129E  00 

3 

0.44370878E 

01 

0.18522095E  00 

4 

0.59226222E 

01 

0. 17510626E  00 

5 

0. 74 130419E 

01 

0. 16357016E  00 

6 

0. 89088589E 

01 

0. 15132422E  00 

7 

0. 10410295E 

02 

0.13896375E  00 

8 

0. 11917337E 

02 

0. 12693599E  00 

9 

0. 13429796E 

02 

0.  11554220E  00 

10 

0. 14947371E 

02 

0. 10495916E  00 

11 

0. 16469693E 

02 

0.  95268123E-01 

12 

0. 17996365E 

02 

0. 86483033E-01 

13 

0. 19526988E 

02 

0. 78574225E-01 

14 

0. 21061176E 

02 

0.71486416E-01 

15 

0.22598569E 

02 

0.  65151455E-01 

16 

0.24138836E 

02 

0. 59496772E-01 

17 

0. 25681676E 

02 

0.  54450669E-01 

18 

0.27226821E 

02 

0. 49945374E-01 

19 

0.28774030E 

02 

0. 45918596E-01 

20 

0.30323089E 

02 

0.42314147E-01 

21 

0.31873811E 

02 

0. 39Q81971E-01 

22 

0.33426028E 

02 

0. 36177855E-01 

23 

0.34979591E 

02 

0. 33562960E-01 

24 

0.36534371E 

02 

0. 31203293E-01 

25 

0.38090251E 

02 

0. 29069161E-01 

26 

0. 39847129E 

02  - 

0.27134649E-01 

27 

0.412  0491PE 

02 

0. 25377148E-01 

28 

0.42763O23E 

02 

0. 23776920E-01 

29 

0.44322885E 

02 

0. 22316715E-01 

30 

0.45882936E 

02 

0.20981445E-01 

LAMBDA  SUM  EQUALS 

0. 24908018E 

01 

TEST  -  0. 45553093, E  02 
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TABLE  #7 


\<BT=1) 

0,  240 102 94 E  00 
0.  644152 16E-01 
0.  35732475E-02 
0.  670133S7E-04 
0.  67770436E-06 
0. 42760217E-08 
0,  18594458E-10 


A.(BT=2) 

0. 30820594E  00 
0. 23550007E  00 
0. 76526055E-01 
0.  77430834E-02 
0. 33492519E-03 
0.  86126762E-05 
0. 15088441E-06 
0.  19268016E-08 
0.  18759706E-10 


Al(BT=3) 

0.31381055E  00 
0. 30429152E  00 
0. 23017193E  00 
0.  82507648E-01 
0. 10966672E-01 
0.  69994684E-03 
0. 284163  34E-04 
0.  82680435E-06 
0. 18234746E-07 
0.  31610705E-09 
0. 44243049E-11 


A.(BT=4) 

0.3 14 1404 IE  00 
0. 31339271E  00 
0. 30140125E  00 
0. 22674445E  00 
0. 86290125E-01 
0. 13513561E-01 
0. 10927087E-02 
0. 58757207E-04 
0. 23453874E-05 
0. 73016896E-07 
0. 18285326E-08 


FLAT  SPECTRUM  EIGENVALUES 
NORMALIZED  SO  THAT  £  =  T 

(SEE  EQUATION  #83) 
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A  common  type  of  digital  communication  system  is  binary  frequency  shift 
keying  (FSK)  whereby  every  T  seconds  the  transmitter  sends  a  pulse  of  one  of 
two  frequencies.  The  receiver  makes  a  decision  (every  T  seconds)  as  to  which 
frequency  was  transmitted.  A  sub  -optimum  receiver  for  this  case  obtains 
estimates  of  the  two  noise  waveforms  by  passing  received  signals  through 
filters  centered  at  the  sending  frequencies  and  then  crosscorrelates  these 
estimates  with  the  received  waveform.  Two  slightly  different  versions  of 
this  cross-correlator  were  considered,  and  the  probability  of  error  for  each 
case  was  calculated.  The  results  seem  to  agree  with  previous  experimental 
work  by  Cossette  and  Wolf. 
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